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In this paper, the current distribution on a bare conducting loop, sit- 
uated in free space over a semi-infinite medium, is obtained for arbitrary 
time harmonic excitations. The loop is assumed to be thin, perfectly con- 
ducting and the standard one-dimensional integral equation and its Fourier 
series solution are used as the starting points. The field due to the current 
in the loop, where the semi-infinite medium is absent, is expressed as a 
superposition of plane waves. The tangential component of the field reflected 
by the interface, of the semi-infinite medium, is evaluated using appropriate 
Fresnel reflection coefficients. This reflected field serves as a new source 
for the loop and induces a current on the loop. The field due to the induced 
current is treated in the same manner, and this process is repeated indef- 
initely. The summation of the original current and all the induced currents 
gives the steady-state current on the loop. 

I. INTRODUCTION' 

It is well known that a high altitude nuclear burst generates an 
intense electromagnetic transient which covers a large geographical 
area. 1 This transient field induces currents in communication circuits 
and, if these are large enough, adversely affects communication channels. 
One problem of particular interest in land-line communication is the 
coupling to large loops formed by cables. 

The loops formed by cables deployed in practical communication 
systems are very complex and cannot be analyzed easily. Typically, 
they run for many miles over inhomogeneous terrain and contain many 
junction points; nevertheless a great deal of insight into the behavior 
of these irregular loops can be obtained by studying the behavior of 
a large regular loop over a homogeneous ground. In this paper, the 
theoretical foundations for an analysis of a circular loop over a homo- 
geneous ground are developed, for time harmonic excitations. The 
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response of the loop for transient fields may be obtained by standard 
Fourier transform techniques. 

Various problems related to the thin circular loop have been con- 
sidered by numerous authors. These may be broadly classified under 
two categories: loop in an infinite homogeneous medium and loop in 
a stratified medium. 

In the first category, Pocklington, Oseen, Hallen, Storer, and Wu 
have analyzed the problem of a bare thin perfectly conducting circular 
loop in free space. 2-6 All these authors use the Fourier series expansion 
to solve the integral equation for the current in the loop. A good analysis 
of the problem is given in Wu's paper. Adachi and Mushiake analyze 
the same problem by solving the integrodifferential equation for the 
current using an iterative method. 7,8 Mei, Baghdasarian and Angelakos, 
and Tang have discussed the direct numerical solution of the integral 
equation. 9 " 11 A variational approach for determining the scattering 
cross section of a loop is given by Kouyoumjian. 12 Problems concerning 
loaded loops are considered by Iizuka, Harrington and Ryerson, and 
Harrington and Mautz. 13-15 The analysis of a loop in a conducting 
medium is an extension of the analysis of a loop in free space and lends 
itself to certain approximations. Kraichman, Chen and King and 
King, and Harrison and Tingley have discussed the bare loop in a dis- 
sipative medium; 10-18 Galejs has discussed an insulated loop in a dis- 
sipative medium. 19 Finally, the solution to the problem of two identical 
coaxial coupled loops in a homogeneous medium has been solved by 
Iizuka, King and Harrison. 20 

In the second category the literature is mostly on small loops or 
magnetic dipoles over different types of media and is very extensive. 
(See Ref. 21 for an extensive bibliography.) Wait has considered the 
problem of loops over a homogeneous earth; 22,23 recently, Sinha and 
Bhattacharya have analyzed the problem of a vertical magnetic dipole 
buried inside the earth. 24 

The treatment of a small current carrying loop as a magnetic dipole, 
while satisfactory for many purposes, is nevertheless inexact. Moreover, 
we do encounter situations where the loop diameter is comparable to 
the wavelength of the excitation frequency and here we cannot assume 
the current to be uniformly distributed on the loop. A typical example 
would be the excitation of a loop by a narrow electromagnetic pulse 
which contains a broad spectrum of frequencies. The purpose of this 
paper is to solve the problem of a bare loop over a homogeneous earth 
taking into account the current distribution on the loop. 

Specifically, the system under consideration is a bare, thin, perfectly 
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conducting circular loop of mean radius b, formed by bending a cy- 
lindrical wire of radius a, situated in free space with its plane parallel 
to and at a distance d from the interface of a semi-infinite, linear homo- 
geneous, isotropic medium (Fig. 1). The loop is excited by an electro- 
magnetic wave of harmonic (exp jut) time variation (the slice generator 
used to compute the admittance being a limiting case). It is assumed 
that k a « 1 and a <<C b, so that, if the loop were situated in a homo- 
geneous medium the current distribution induced by a specified time 
harmonic excitation is given by the so-called one-dimensional integral 
equation. In addition it is assumed that d»fl. 

It is desired to determine the current distribution /(<£) on the loop 
in the aforementioned system. This is accomplished in the following 
three, more or less self-contained, sections. In Section II the field of a 
circular filamentary current in free space is expressed as a superposition 
of plane waves. Section III evaluates the reflected field when an arbitrary 
field, of the general form obtained in Section II, is incident on the 
interface between free space and the semi-infinite medium. In Section 
IV the results of Sections II and III are combined with the integral 
equation for the current on a loop in free space to determine the steady 
state current 1(4) on the loop using a recurrent "reflection-induction" 
scheme. 



II. THE ELECTROMAGNETIC FIELD OF A CIRCULAR FILAMENTARY CURRENT 

Consider a circular filament of current I(4>) (Fig. 2a) of radius b 
situated in free space. The coordinate system is so chosen that the 
loop is parallel to the xy plane at a distance d from it. The loop current 
may be expressed as a surface current density K, in the z = d plane, 
in the following manner. 

K = a#/fo)5( P - b)5(z - d) (1) 

where p, <f> and z are the cylindrical coordinates and a^ the unit vector 
in the <f> direction. The electromagnetic field due to K may be expressed 
as a superposition of plane waves as follows. 25 

H ( ° = ff* [±P, ±Q, (IP + mQ)(l - I 2 - m 2 )- § ] 

•exp [jk [lx + my =F (1 - I 2 - m 2 ) h (z - d)]\ dl dm, 



z <t d. (2) 
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Fig. 1 — A thin circular loop over a homogeneous semi-infinite medium. 

E< ° = v ° \Jl {[lmP + (1 " l2)m " f " m2) "' 

- [(1 - m 2 )P + MHU - Z 2 - m 2 )"*, ±(IQ - mP)} 
•exp {jfc [te + my =F (1 - Z 2 - m 2 )*(z - d)]) dZ dm, 



2 ^ d, (3) 
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Fig. 2 — (a) A filamentary loop current in free space, (b) A filamentary loop current 
over a homogeneous semi-infinite medium. 
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where 
P(l, m) - (k 2 /8r 2 ) JJ ^ K u (x, y) exp [-jk (lx + my)] dx dy, (4) 

Q(l, m) = -(k 2 /87r 2 ) JJ_ K x (x, y) exp [-jk (lx + my)] dx dy, (5) 

k 2 = co 2 fi e and -q = (n /e ) h . (6) 

The quantities I and m are in general complex and the integrals in 
equations (2) and (3) are contour integrals, the choice of contours 
being dictated by physical considerations. A possible choice of I and m 
is defined by the following transformation 

I = T cos \I/,\ n _ , , . 

\ ^ t, y being complex. (7) 

m = t sin \f/,\ 
Let 

+ 00 

/(*) = E I n exp(jn<j>). (8) 

Substitution of equations (1), (7) and (8) into equations (4) and (5), 
and changing from rectangular to cylindrical coordinates yields the 
following equations 

P(r, rf) = (klb/Sw) £ (- j) n exp (jn^J^kMKh^ + /„+,), (9) 

T1--00 

+ 00 

Q(t,*)= -Kklb/Sir) 2 (-jTexp(Jnf)J n (kM)(.In-i-In + i), (10) 

n = — co 

where J H denotes the Bessel function of order n. 

III. CIRCULAR FILAMENTARY CURRENT OVER A SEMI-INFINITE MEDIUM 

Here again we consider the filamentary loop of Section II, but instead 
of being situated in free space it is situated over the homogeneous 
semi-infinite medium z < (Fig. 2b). The total field in the region 
z^O consists of the primary field E (,) , H ( " of the filamentary current 
[equations (2) and (3)] and the field E (r> , H (r) reflected by the interface 
2 = 0. We proceed as follows to evaluate the latter. Let 

H (,) (.-c, y, z) = J J H a (l, m) exp [-jk^il, m)-r] dl dm, z ^ d, 

(11) 
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E (0 (x, y, z) = Jf E„(Z, m) exp [-jkji (l, w)t] dZ dm, z ^ d, 

(12) 
where 
H„(J, m) = [-a,P - a„Q + a,(/P + mQ)(l - / 2 - m 2 ) _i ] 

•exp[-iM(l - I 2 ~ m 2 f], (13) 

8.(1, m) = ij {a x [ZmP + (1 - Z 2 )Q](1 - f - m 2 )"* 

+ a„[(m 2 - 1)P - JmQ](l - I 2 - m 2 )" 5 + a.(mJ° - IQ)} 
■exp[-jk d(l - I 2 - m 8 )*], (14) 

n (/, m) = -aj - a„ra - a z (l - I 2 - m 2 ) h , (15) 

and 

r = a,x + a„?/ + a r z. (16) 

Each one of the constituent plane waves propagates in the "direction" 
n.(l, m). Let R x (l, m) and RS> m ) represent the Fresnel reflection 
coefficients for the cases where the incident electric field is perpendicular 
and parallel respectively to the plane of incidence. Evaluation of the 
reflected field is achieved by resolving each plane wave into components 
with the electric field perpendicular and parallel to the plane of incidence. 
To this end we define a local coordinate system 1 as shown in Fig. 3. 

The plane of incidence is defined by the unit vectors a, and n„(Z, m). 
Let 

a, = (a, Xn )/[1 - (a z -n ) 2 ] 5 , (1?) 

a2 == a* }\ aj . 

Then a, , a x and a 3 form a right-handed coordinate system, &i is normal 
to the plane of incidence and a, and a 2 lie in the plane of incidence. 
In terms of a x and a„ we have 

&l = (I 2 + mVW - &„), (18) 

a 2 = (I 2 + m 2 )-'(/a x + ma,). (19) 



t The propagation vector n becomes complex for certain values of I and m, the 
associated plane wave being inhomogenous. When this happens some of the terms 
used in the analysis (for example, plane of incidence, coordinate system, normal, 
and so on) become inaccurate and should be interpreted in a generalized sense. 
The results obtained are quite general and valid. 
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PLANE OF INCIDENCE 
CONTAINS — — 

n«,8 



INCIDENT WAVE 
DIRECTION fi ) 




CCLj PLANE (INTERFACE) 
CONTAINS 



Fig. 3 — Local coordinate system. 

The field vectors H„ and E„ are now resolved into components per- 
pendicular (H 0j . , E.j.) and parallel (H oB , E „) to the plane of incidence 
as follows 



H„ — H oX -j- H || , 



(20) 



where 



(22) 
(23) 



H oX = a,(Z 2 + m*)~\lQ - mP) exp [-jk d(l - I 2 - ro 2 )*], (21) 
H „ = [a,(l - I 2 - m 2 )"* - a 2 (l 2 + mY*](lP + mQ) 
■exp[-jk d(l - I 2 - m 2 )*], 
E = E oX + E oll , 
where 
E oX = •,*.(? + m 2 y\l - I 2 - mTKlP + mQ) 

exp [-jk d(l - I 2 - m 2 ) 1 ], (24) 
E„„ = ,.[a, - a 2 (/ 2 + m 2 )-»(l - I 2 - m 2 )%nP - IQ) 

•exp [-jk„d(l - I 2 - m 2 )*]. (25) 
It may be easily verified that 

H oU = v : l (n„ X E.J.), (26) 

E„ = -,„(n. X H.J. (27) 
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The reflected field corresponding to the incident field defined in 
equations (11) and (12) may be represented as 

H (r) (z, y, z) = fj" H 2 (Z, m) exp [-jkMh m)-r] dl dm, z ^ 0, 

(28) 

E <r> (a;, y, z) = \\ E 2 (l, m) exp [-jkji 2 (l, m)-T] dl dm, z ^ 0, 

(29) 

where 

n 2 (Z, m) = -aj - a„m + (1 - I 2 - w a )*a. . (30) 

Let 

H 2 = H 2 ± + H 2 |) , (31) 

E 2 = E 2X + E al . (32) 

According to Fresnel's laws 

E 2X = RS, ™)E oX , (33) 

H 2X = R t (l, m)H oX , (34) 

where 

R L (l, m) - 1/Jb.d - I 2 - m 2 )* - M [fc 2 - k\{l % + ^ 2 )] i l/ 

{ M fc (l - I 2 - m*)* + n [k 2 - kl(l 2 + m 2 )]'}, (35) 

R t (l, m) = { Mo fc 2 (l - I 2 - m 2 )* - fxk [k 2 - kl(l 2 + m 2 )]*}/ 

{ Mo fc 2 (l - I 2 - m 2 f + nk [k 2 - kl(l 2 + m 2 )]*}, (36) 

A; 2 = —join(a- + Jttc). (37) 

We also have the relations 

H 2I1 = ^'(n, X E 2X ), (38) 

E 2U = -rj.fa X H 2X ). (39) 

Substitution of equations (33), (34), (38) and (39) into equations (31) 
and (32) yields 

E 2 = ft x E oX - i7„ft,i(n 2 X H oX ), (40) 

H 2 = R t K oi . + v^Rifa X E oX ). (41) 
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Equations (40), (41), (28) and (29) specify the reflected field com- 
pletely. The 4> component of the electric field is of special interest since 
it induces a current in an actual circular loop. Taking the <p component 
of equation (40), using polar coordinates for x, y, and equation (7), 
we obtain 

#2* = i?.[B|(t)(1 - r 2 )*(Psin $ - Q cos $ sin (</- - tf>) 

- R ± (t)(1 - t 2 )~ 1 (P cos t + Q sin $) cos (^ - <*>)] 

•exp[-;Av/(l - r 2 ) 5 ], (42) 

where 

B,(t) = W(l " r 2 ) } - ,*.<*» - /C 2 r 2 ) 4 ]/ 

W(l - r 2 )* + nK(k 2 - klr 2 f], (43) 

B x (t) = [ M A-„(1 - t 2 ) 5 - vAk* - A: 2 „r 2 )']/ 

[ M /c (l - r 2 ) 5 + Mo (fc 2 - fc'oT 3 )*]. (44) 

P(t, ^), Q(t, yf/) are defined in equations (9) and (10), and after re- 
arrangement yield the following equations 

P cos $ + Q sin ^ = i(47r) _1 A: 2 & £ U'JJcM) exp jn(^ - |) , 

(45) 

Psin ^ — Q cos ^ = (4TT) _1 n/c„ ^ I n Jn(k bT) exp j»(^ — -J 

n--co V -/ 

(46) 

where the prime denotes differentiation with respect to the argument. 
The component of the reflected field is contained in equation (29) 
and is explicitly given by 



E l ;\ P , <f>, z) = f f E 2 ,(r, +) 
■'o J c 



•exp [jk n pT cos (\J/ — 4>) — jk n z(l — t 2 )'] d\f/r dr. (47) 

The contour of integration C in the complex \f/ plane is to be chosen 
from physical considerations. An examination of equations (42), (45) 
and (46) shows that the \f/ integrals gi , <0 2 are of the form 



<0i .2 - j c ™ (<£ ~ 0) exp \jn\f - |J + ;'A- n pr cos (^ - 0) J drf/. 



(48) 
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Since $ lti are solutions of Helmholtz equation in cylindrical coordinate 
system, we expect them to be cylindrical functions (compare with Ref. 
26, pp. 367-368). The requirement that they be bounded when the 
argument of the Bessel functions approach zero determines that they 
are the ordinary J functions. Thus we obtain, 

& = 27m(fc pT)"V n (fc pT) exp (jn<p), (49) 

& = -j2irJ' n {k oP T) exp (jn<t>). (50) 

Substitution of equations (42), (45), (46), (48), (49) and (50) into 
equation (47) yields the following expression for the component of 
the reflected field 

+ 00 

K\p, *, *) = 2tt E IMp, z) exp (jrup), (51) 

T1--O0 

where 
z n ( P ,z) = (4x)- 1 /c 2 &r ?0 P [n\k 2 oP b)- l J n {k oP T)J n {kM)r- l iX - t')%(t) 

• 

- J' n {k bT)J'n{k oP T)r{\ - T 2 y h R x (r)] 
•exp [-jk (z + d)(l - r 2 ) J ] dr. (52) 

In particular, 

z n (b, d) = (feT'ftfrr. r [7i 2 (fc„6)-V*(*; 6T)T- l (l - 0%(r) 

- J' n \k br)r{l - 0"*B x (t)] exp [-J2M(1 - t 2 ) 1 ] dr. 

(53) 

IV. CONDUCTING CIRCULAR LOOP OVER A SEMI-INFINITE MEDIUM 

In this section we consider a perfectly conducting thin circular loop 
of mean radius b situated over the semi-infinite medium (Fig. 1) with 
its plane parallel to the interface. The radius of the wire forming the 
loop is a and the height of the loop above the interface is d. 

4.1 Current Distribution 

Let El 0) (<p) be the <p component of the applied electric field, I (o) (<p) 
the current distribution that would be created by J2{ 0) (<p) on the loop 
if it were situated in free space and 7(0) the current distribution caused 
by Ei* (0) when the loop is situated as shown in Fig. 1. Typical examples 
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of El o) (<t>) are the slice generator (used in admittance computations) 
and the <£ component of the electric field that would exist at the loop 
location in the absence of the loop (as in scattering problems). 

The relation between the applied tangential electric field and the 
currents induced on a loop takes the form of coupled integral equations 
which are extremely difficult to solve. 6 However if the loop is "thin" 
(a « b, a « X) the current distributions on the loop is given accurately 
by the so-called "one-dimensional integral equation." Thus we have 

El°\4>) = r G(<f> - *')7 (0 V) <fc', (54) 

where 
G(x) = X4T)- 1 v.\k.b cos x + (Kb)- 1 gj W 

• f \R(X) eXp [- jkMx)] }a^2a .in e dd , (55) 

R(x) = 6[4sin 2 (x/2) + (a 2 /b 2 )]K (56) 

A formal solution to equation (54) is obtained by using Fourier series 
representations as follows. Let 

/ ( " > (0) = 2 I? exp(^), (57) 

K\4>) = 2 «n (0) exp (jruj>), (58) 

71= -SO 

G(<f> -<*>') = £ Pn exp \jn(<l> - <*>')], (59) 

n= — w 

where 

I { n 0) = (27T)" 1 f T /(<*>) exp (-jn<f>) d<t>, (60) 

■'0 

a ( n 0) = (27T)" 1 f ' El'\4>) exp (-jn<t>) <ty, (61) 

•'0 

ft = (2t) _1 [ * G(x) exp (- jnx) dx. (62) 

Substitution of equation (57), (58) and (59) into equation (54) yields 
the following relation between the Fourier coefficients 

Tl n) = (2 7 r)- , (al 7/3 n ). (63) 
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The current distribution is given by 

J (D, (4>) - (27T)" 1 E U 0) M exp(in0). (64) 

n = — oo 

Let E+* (0) be the tangential electric field, at the loop, of the reflected 
field whose incident field is caused by i" (o> (0) and let I 0) (<f>) be the 
current distribution on the loop caused by #"'(<£). Let 

E? (*) = E «i 1} exp (**), (65) 

/ (,, (*) - E Jj" «*(**). (GG) 

n=-=o 

Then 

J< 1) = (2 T )- 1 [ai 1> //3»] (67) 

[compare with equation (63)]. Also from equations (51) and (65) we 

obtain 

q£» = 2*I ( n 0) z n (b, d). (68) 

Hence 

/" } = In'Mb, d)//3j. (69) 

In general, if Zj M denotes the ?ith Fourier coefficient of current distribu- 
tion 7 <A:) (</>), induced by the fcth reflected field we have 

/<*> = /'^"[^(b, cO/jBj, /c = 1, 2, • • • . (70) 

Let 

/(*) = E /nexp(^). (71) 

n- — oo 

Then 

/„ = E/r, (72) 

that is, 

in = n o) + i? +/? + ••■. 

= I^ { 1 + [«.(&, d)M + [*»(&, fl/ftO* +•■•}, 

= Jj #, U - [*»(&, (DM)' 1 , provided | [z n (b, d)/AJ \ < 1. (73) 
Henceforth, we assume that | z n (b, d)/p n \ < 1. Substituting for I m mi 
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in equation (73) from equation (63) we obtain 

/„ = (27rn*n/3n -*»(&, <or\ (74) 

/(<*>) = (IM- E aj 0) [ft, - z n (b, d)]- 1 ex V jn<t>. (75) 

n=> — co 

The following expression for /3„ is obtained by simplifying equation 
(62): 

ft. = i(47r/c„6 2 )- , r 7 „(2/7r) 

• P (*(*.&)'[#.«(«) + -l/ n -,(a)J - 7i 2 M n (a)} fl -. 2aain9 </0, (76) 

Jo 

where 

M B (x) = - f 2 (sin 2 + .r74bT 5 

•exp [-i2/c 6(sin 2 6 + x 2 /4b 2 ) h ] cos (2r?0) dd. (77) 



Let 

r/2 



7T. = (2/r) f 2 .l/„(a) 

->o 



dd. (78) 

o-»2a sin 9 



Then 

ft, = i(47r/c 6 2 )- 1 77o [KA- &) 2 (/v n+1 + #„_,) - n"JKJ. (79) 

■ 

4.2 Inpuf Admittance 

Let the primary source be a slice generator (delta function source) 
of voltage V located at <t> = 0. That is 

E ( ;H<f>) = [V5(<t,)/b]. Gso) 

Substituting equation (SO) into equation (61), we obtain 

„<-' = (7/2x6). (81) 

Substitution of equation (81) into equation (75) yields 

I® = 7(4/6)"' j£ [ft. - z n (b, d)}- 1 exp (jn4>). (82) 

The admittance Y at the input terminals at <f> = is given by 

F = I(0)/V = (4tt 2 6)- 1 £ [ft, - z n (b, d)]- 1 . (83) 
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The use of the delta function generator will give rise to an infinite 
input admittance 8 so that the series in equation (83) is divergent. 
However this difficulty is overcome by computing the difference be- 
tween the admittances of two loops of different radii. 27 

V. SPECIAL CASES 

5.1 The Magnetic Dipole Over a Semi-Infinite Medium 

When the radius, b, of the loop becomes very small compared to the 
wavelength (that is, k b « 1) the current distribution on the loop 
becomes uniform. This enables us to retain only the zeroth terms in 
the infinite series representing the different quantities of interest. The 
field of a dipole over ground is well discussed in the literature and will 
not be considered here. The input impedance of a dipole over a semi- 
infinite medium is of considerable interest and may be obtained from 
equation (83). Thus we obtain 

Z ia = 47r 2 &[/3„ - z (b, d)]. (84) 

The term 47r 2 5/3 represents the input impedance of the loop in the 
absence of the semi-infinite medium and the term — 4:ir 2 bz represents 
the contribution of the semi-infinite medium. That is 

Z ia = Z pri + Z BO0 , (85) 

where 

Z Dt , = 4/6/3, , (86) 

Z aeB = -4**bz.(b, d) 

= T(k b) 2 ri a f M J'i(kobr) T (l - T 2 y h Rx(r) 
Jo 

•exp[-j2k d(l - T) h ]d T . (87) 

5.2 Thin Circular Loop Over a Perfectly Conducting Plane 
Let 

z*(jb,d) = limz n (b,d). (88) 

1T-«00 

When a — > » the reflection coefficients simplify to 

i2,( T ) = +1, R ± (r) = -1- (89) 

Substituting equation (89) into equation (53), we obtain 
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Z*(b, d) = (47r)- 1 ^6r 7o 

• f" [n 2 (k by 2 J 2 n (k e bT)T~\l - r 2 )* + J^ 2 (fc 6r)r(l - t 2 )"*] 

.exp[-i2fc d(l - T 2 ) l ]dT. (90) 

The integral in equation (90) may be simplified, by making use of the 
properties of Bessel functions, and yields 

z*(b, d) = (47r&)- l U£(W>) 2 (f n -i + f-i) - fftj, (91) 

where 

I". - f *G - r 2 )"V 2 (/c 6r) exp [-i2M(l - t 2 )*] dr. (92) 

«>o 

An alternate expression for f„ is obtained as follows: 

j2(fc o 6 T ) = - f J (2k b T sm 6) cos (2nd) dd; (93) 

substituting equation (93) into equation (92) and changing the order of 
integration we obtain 

Tn = j(k n by 1 M„(2d) > (94) 

where M n (x) is defined by equation (77). Therefore, 

Z*«Q>, d) = i(47rfc 6 2 )- 1 

•Vo\(h)(k b) 2 [M n+l (2d) + M n - l (2d)] - n 2 M n (2d)\. (95) 

A comparison of equations (95) and (76) reveals a strong similarity 
between the expressions for /3„ and z* n (b, d). The input admittance of 
a thin circular loop over a perfect ground plane is given by 

Y = (4/6)"' £ ft - z*(b, d)]-\ (96) 

The above formula agrees, with that derived by Iizuka and others, 20 
for the input admittance F (a> of a loop in the presence of an identical 
coaxial loop carrying a current distribution which has an opposite 
phase. They use the simpler Kernel given by Storer 5 to computer /3 n , 
a procedure satisfactory for small loops. However they make use of the 
similarity between /3„ and z*(b, d) to compute the latter using the 
approximate expressions given by Storer. This procedure will yield 
erroneous results for large separation, d, for the following reason. 

The approximate expressions for /3„ given by Storer, or for that 
matter Wu, are valid for k a a <3C 1. The corresponding condition to be 
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imposed in the evaluation of z*(b, d) is, (2k e d) « 1. Thus it is seen that 
the approximate expressions of Storer give accurate values of z*(b, d) 
only for very small separations. 

VI. NUMERICAL COMPUTATIONS 

Equation (76) was further analyzed and approximate expressions 
for the /3„ coefficients were derived in terms of Bessel and Legendre 
functions. The integral defining the z n coefficients could not be expressed 
in terms of known functions and so was evaluated by numerical integra- 
tion. However, in the evaluation of z* it was possible to use some of 
the formulae developed for /3„ for small values of d [compare with 
equation (95)]. The numerical integration was carried out by using 
the Romberg integration scheme. All the computations were done by 
Fortran programs on a GE-635 computer. 

Figure 4 shows the variation of the input admittance of a loop, 
when it is in free space, over moist earth and over an infinitely con- 




0.8 0.9 1.0 

(2^rb/X) 

Fig. 4 — Input admittance of a circular loop, over different media, as a function of 
frequency (d/b = 0.25; a/b = 0.002). G— conductance; B— susceptance. I— loop in 
free space; II — loop over moist earth; III — loop over perfect ground (<r — * <*> ). 
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ductive ground plane, as a function of frequency. The values of the 
various parameters used in these calculations are 2irb = 30 meters, 
/ = 5 MHz to 13 MHz, d/b = 0.25, a/26 = 0.001, <r = 5 millimhos/ 
meter, and e/e = 15. The last two parameters characterize the moist 
earth. The frequency range was deliberately chosen so that the moist 
earth cannot be approximated either as a highly conductive medium (low 
frequency approximation) or as a lossless dielectric (high frequency 
approximation) . 

The real part G of the input admittance shows its characteristic peak 
near k b = 1 (these occur for values of k b near 1, 2, 3, • • •) but the 
exact location of the peak as well as its magnitude depends on the 
medium below the loop. When the loop is located above a highly con- 
ducting ground plane the resonance is particularly sharp at k b = 1 
since at this frequency the loop and its image are exactly half wavelength 
apart. The imaginary part B of the input admittance changes from 
inductive to capacitive near k b = 0.7 and back to inductive near 
k b = 1. Here again the transition at k b = 1 for the loop over a highly 
conducting ground plane is almost discontinuous. 

Figure 5 shows the variation of the input admittance of a loop over 
a highly conductive ground plane as a function of the distance between 
the loop and the ground plane. The curves are plotted for k b = 1, 
a/26 = 0.001 and d/b ranging from 0.25 to 5.0. It is observed that as d/b 
increases, the input admittance approaches the free space value in an 
oscillatory manner. 

The aforementioned calculations are presented only as examples 
of the different types of investigations that may be carried out based 
on the theory developed. The computer programs developed in this 
connection are very general and may be used for computations of loops 
as large as k b = 10. 

VII. SUMMARY 

The problem of a thin, perfect]}- conducting, circular loop situated 
in free space over a semi-infinite homogeneous isotropic medium was 
solved. Expressions for the current distribution on the loop caused by 
an arbitrary time harmonic source [equation (75)] and the input admit- 
tance [equation (83)] were derived. The results are applied to special 
cases to evaluate the input impedance of a vertical magnetic dipole 
over a semi-infinite medium [equation (84)] and the input admittance 
of a circular loop over a perfectly conducting ground plane [equation 
(96)]. Some numerical results are also given. The analysis for the general 
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Fig. 5— Input admittance of a circular loop over an infinitely conducting ground 
plane as a function of height (*o& = 1, a/b = 0.002). G— conductance; B— suscept- 



case of a circular loop in an arbitrary homogeneous medium (as opposed 
to free space assumed here) over another homogeneous medium can 
be easily done by modifying the parameters. 
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